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ABSTRACT 

We present the results of spectral analysis of 42 simultaneous broad-band Ginga-OSSE and 
RXTE-OSSE observations of Cyg X-l carried out in 1991 and 1996-1999. The broad-band 
spectra from 3 to ~ 1000 keV can be well described by thermal Comptonization model with 
reflection from the cold disc, with an additional soft component visible below 10 keV. The 
relative contribution of this component to the total energy flux appears to be higher in the 
spectra with larger reflection amplitude and steeper photon index of the thermal Comptonized 
component. We consider a number of physically realistic models to describe the shape of the 
E<10 keV excess. The additional soft component can result from thermal Comptonization 
by electrons with a low Compton parameter, or can be a part of a nonthermal, power-law like 
emission extending above 1 MeV. 

We study correlations between parameters obtained from the spectral fits with different 
models. We confirm a general correlation between the photon index T and the amplitude of 
reflection R. We find that simple phenomenological models (like power-law plus Compton 
reflection) applied to the narrow band (3-20 keV) data overestimated the values of R and V, 
although the simple models did rank correctly the spectra according to R and T, as it was 
demonstrated in the original publications on this subject. 

The dynamic corona model provides a satisfactory description of the observed correla- 
tion, while the hot inner disc models have problems in reproducing it quantitatively. On the 
other hand, in the context of the dynamic corona model it is difficult to understand correlations 
with the timing characteristics, which seems natural in the hot disc scenario. We do not find 
significant correlation between the electron temperature and other spectral parameters, while 
the optical depth of the hot medium seems to decrease when the spectrum becomes softer. It 
is also shown that spectral parameters are well correlated with the timing characteristics of 
the source. 

Key words: accretion, accretion discs - black hole physics - gamma-rays: observations - 
stars: individual: Cygnus X-l - X-rays: binaries 



1 INTRODUCTION 

Matter accreting onto a black hole, whether supermassive in Seyfert 
galaxies or stellar mass in Galactic X-ray binaries, releases most 
of its gravitational energy in the form of X-rays deep in the po- 
tential well. Accretion may proceed in a number of regimes. The 
accreting gas can approach the black hole in a disc-like con- 
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figuration fchakura & Sunvaevlll973l) if the gravitational energy 
is effectively transported away in the form of radiation, or in a 
form of an almost spherical flow if the energy exchange mech- 
anism between pro tons (carrying most of the energy) and elec- 
trons is inefficient ISha piro. Lightman & Eardleviri976t Echimaru 
Il977t l Naravan, Mahadevan & Ouataert 1998). Magnetic fields can 
play an important role transporting large fraction of the total avail- 
able energy and dissipating it in a rarefied medium (corona) above 
the disc iGaleev. Rosner & Vaianal Il979t iTout & Pringld Il992t 
ISvensson & Zdziarskilll994 lBeloborodovlll999d: iMiller & Stonel 
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2000). However, every model is based on a number of assumptions 
that include prescriptions for the viscosity, the vertical distribution 
of the energy release through the flow, the energy transport mecha- 
nisms, etc. Given the difficulties in the accretion physics, observa- 
tions should help in choosing among the different possibilities as 
well as guide theoreticians in the right direction. 

Cygnus X-l, one of the best studied black hole binaries 
(BHB), served as an accretion disc laboratory since the end of the 
1960's. The most dramatic observed phenomena are the spectral 
state transitions occurring every few years, when the source, typ- 
ically emitting most of its energy at about 100 keV in the hard 
state, switches to a soft state consisting of a prominent ~ 1 keV 
black-body and a power-law-like tail. The hard state spectrum was 
believed to orig inate from thermal Co mptonizati on in a hot elec- 
tron cloud IShap iro. Lightman & Eardlevl I"l97(& llc himarul fT977t 



ISunvaev & Triimpeill979tlSunvaev & Titarchukll980h . The black- 
body looking soft-state spectrum was associated with the optically 
thick accretion disc l Shakura & Sunva evl 19731) . while the origin of 
~ 100 keV emission in that state was not discussed much, because 
the detailed spectrum was not available. 

During the last decade, the quality of the XAy-spectra in- 
creased dramatically improving our knowledge as well as pro- 
ducing many new questions. In the hard state, the spectrum turns 
out to be rather complicated containing a number of compo- 
nents. Thanks to the broad-band coverage by Ginga and CGRO 
(and later by ASCA, RXTE and BeppoSAX), and advances in 
modell ing of Com ptonizat ion at mildly r elativistic temperatures 
lCoppilll99l IPoutanen & Svenssoni ri996). the parameters of the 
electron cloud where Comptonization takes place were deter- 
mined to a high accuracy. In Cyg X-l, the electron tempera- 
ture of kT c ~ 100 ± 50 keV and Thomson optical depth of 
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flection bump, a signature of the presence of cold matter in the 
vicinity of the X-ray emit t ing source, was discov ered I Do ne et alj 
ll992tlEbisawa et alll99^ : lGierlinski et alll997l). The black body 
associated with the cooler accretion disc (Balucinska & Hasinger 
ll99ltlBamcihska-Church et alJll995tlEbisawa et alJll996h and an 
additional soft excess at a few keV of unknown origin furthe r com- 
plicat e spectral decomposition iDi Salvo et alj2001rlFrontera et alj 
2001). A high energy excess at > 500 k eV discovered by CGRO 
iMcConnell etaflll994lLing et alJI 19971) gives some clues about 
presence of nonthermal particles in the source. 

The XAy-ray soft state spectrum has been studied exten- 
sively by simultaneous observations with ASCA, RXTE, Beppo- 
SAX, and CGRO during summer 1996. In addition to the dom- 
inating black-body, a long power-l aw like tail e xtending up 
to 10 MeV was discovered iMcConnell et alJl2002T) . The high- 
energy spectrum could be well described by single Compton 
scattering off electrons having a nearly power-law distribution 
IPoutanen & Coppilll99^: |Poutanedll998t iGierlinski et aljfT99sl 
hereafter G99; Front era et alj|200lh . The Compton reflection was 
stronger than in the hard state, which was interpreted as a change 
in the geometry of the system, from the hot inner flow in the 
hard state to the standard Shakura-Sunyae v disc with a no nther- 
mal corona in the soft state I Bisnovatvi-Kogan & Blinnikov 1 19771: 
l lchimarul 19771 IPoutanen. Krolik & Rvdell997l : lLi & MilleJl997t 
lEsin et alJll998l : IPout anen & Coppi 1998). Evaporation and con- 
densation of the gas can provide a physical basis for the change 
of th e transition radius between standard disc and hot inner 
flow iMever. Liu & Mever-HofmeistedEoOd : iRozanska & Czernvl 



2000). A smaller reflection fraction in the hard state, how- 
ever, can be explained by a beaming of the primary emis- 
sion away from the disc due to mildly relativistic motion of 
the emitting plasma (ejection model, s ee iBeloborodovl Il999a|pt 
Malzac. Beloborodov & Poutanen 2001, hereafter MBP01). In or- 
der to distinguish among different possibilities, one needs to com- 
pare model predictions with other observational facts. 

A few, well separated in time, broad-band spectra do not 
give us a good picture about relations between different compo- 
nents such as, for example, Comptonized continuum and the re- 
flection bump. One can stud y these relations on a larger data se t 
in a narrower energy band. Zdziarski, Lubinski & Smith 1 1999) 
(hereafter ZLS99) and iGilfanov. Churazov & Revnivtsevl 1 1999) 
(hereafter GCR99) analyzing data from Ginga and RXTE, re- 
spectively, showed that the photon spectral slope of the Comp- 
tonized emission, T, is strongly correlated with the ampli- 
tude of Compton reflection component, R — fl/2n, where 
SI is the solid angle the cool material covers as viewed from 
the source of primary X-rays. This correlation exists for in- 
dividual BHBs and Se yfert galaxies as well as in a sam- 
ple of sources (see a lso Gilfanov, Churazov & Revnivtsey]|200pJ; 
iRevnivtsev. Gilfanov & Churazov! 120011) . IZdziarski et alj 1 2003) 
studied possible statistical and systematic effects and concluded 
that the correlation exists beyond any reasonable doubts. Similar 
correlation also exists for Fourier-frequency resolved spectra, i.e., 
those corresponding to the variability in a given range of Fourier 
frequencies iRevnivtsev. Gilfanov & Churazov 1999). 

The observed correlation provides extremely important clues 
to the geometry of the accreting material and can be used for testing 
theoretical model. The fact that T and R are correlated is a natural 
consequence of co-existence of the cold media (accretion disc) and 
a hot Comptonizing gas in the vicinity of the black hole. The cold 
material acts as a source of seed photons for Comptonization and, at 
the same time, reflects and reprocesses the hard radiation produced 
in the hot gas. 

The cold disc with the hot inner flow model naturally produces 
the cor relation i f ther e is an overlap between hot and cold phases 
IPout anen et alll997l: ZLS99). However, the F-R dependence ob- 
served in BHBs can be quantitatively reproduced only if the ratio 
of the seed photon tempera ture, kT Bee d, to th e electron temperature 
is about 10~ 4 fZLS99: lGiifanov et al.l2 000). For kT c ~ 100 keV 
this gives fcT sco( j ~ 10 eV which is an order of magnitude smaller 
than the disc temperature in BHBs and closer to that expected from 
Seyferts. For kT scc d ~ 300 eV, the spectra are too hard for the 
given reflection fraction. In this model, the spectral slope is an ex- 
tremely steep function of the overlap between t he corona a nd the 
disc, while the reflection varies very little lBeloborodovll200ll) . In- 
trinsic dissipation in the disc can make spectra softer for the given 
reflection, but then the spectral slope will be even steeper func- 
tion of the overlap. Values of reflection larger than 1 sometimes 
observed in Seyferts also cannot be explained. On the other hand, 
all the data ca n be well described by the ejection model with the 
correct fcT aee d tB eloborodovl 1999Jbl: MBP01). 

The problem is that theoretical models try to reproduce the 
best-fitting F and R which are subject to a number of system- 
atic effects. The photon index and, especially, the strength of the 
Compton reflection depend crucially on the spectral shape of the 
underlying continuum and descripti on of the re flection physics 
(see detailed discussion in IZdziarski et ail [2003 ). All the papers 
above, where this correlation was studies, assumed the underly- 
ing spectrum to be a power-law. Since the Comptonization spec- 
trum has a cut-off at high energies it gives fewer incident pho- 
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tons that are available for reflection jWgave r. Krolik & Piedfl998l : 
IPerola et alJl2002HMalzac & Petruccl2002t) . and thus the fitted R 
would be larger. Fitting an exponentially cut-off power-law to the 
broad-band data (e.g. Matt 2001; Per ola et alJl200 2") does not im- 
prove the situation, since this model do es not reproduce well the 
shape of the Comptonization continuum (Malzac & Petrucci 2002; 
Zdziarski et al ] l200l) . Approximate treatment of ionization could 
be another source of errors. 

As we discussed above, the spectra of BHBs are rather com- 
plex having number of components, and it is not possible to re- 
solve different spectral components (e.g. thermal Comptonization 
and soft excess) in the narrow energy range of an instrument such 
as RXTEfPCA or GingafLAC. There is a danger that different com- 
ponents overlap in that energy band producing effectively a power- 
law of one index while in reality the slope of the primary emission 
could be rather different. The resulting reflection amplitude could 
be also influenced significantly. Thus, in order to obtain actual F 
and R to be used in theoretical models, analysis of broad-band data 
with physical models (such as Comptonization) is absolutely nec- 
essary. 

Additional sources of information are the width of the fluores- 
cent Fe Ka line at 6.4 keV and the frequencies of the quasi-periodic 
oscillations that were also obs erved to correlate wit h the reflection 
fraction and the spectral slope iGilfan ov etall200d GCR99). This 
seems to be consistent with the variations of the inner cold disc ra- 
dius. Not much data exist on the variability of the electron tempera- 
ture and Thomson optical depth of the Comptonizing source which 
can provide information about the nature of the emitting plasma 
(electron-proton or electron-positron). It would be of interest to de- 
termine how the optical depth changes with the bolometric flux, 
because this can help in distinguishing the accretion mode the flow 
is in. 

In this paper, we analyze a large set of simultaneous broad- 
band spectra. Four observations of Cyg X-l by GingafLAC and 
CGRO/OSSE from 1991 as well as 38 observations by RXTE/PCA, 
RXTEfHEXTE and CGRO/OSSE from 1996 - 1999 are studied 
in details. For the spectral analysis, we use physically motivated 
Comptonization models and study correlations between model pa- 
rameters such as spectral slope of the primary Comptonization con- 



tinuum, reflection amplitude, the width of the Fe line, electron tem- 
perature of the hot gas, and its Thomson optical depth. 



2 OBSERVATIONS AND DATA ANALYSIS 

The observation log is presented in Table Q Data reduction for 
RXTEhiLS been carried out using LHEASOFT 5.3.1 software; PCA 
responses were generated using pcarsp v. 10.1 and HEXTE re- 
sponses were used from 2000 May 26. In PCA data reduction, all 
5 PCUs were taken into account, when possible. If not all PCUs 
were turned on, we use PCUs 0, 2 and 3. Judging from the Crab 
data, these two PCU configurations produce similar spectral slopes. 
Systematic errors of 0.5% were added in quadrature to the PCA 
data. CGRO/OSSE spectra were prepared by adding per-orbit data 
(average exposure 2-5 ksec), with total exposure up to 12 hours 
and contemporary PCA observation in the middle of the period. 
The systematic error in OSSE spectra varies from 3% at 50 keV 
to 0.3% at 300 keV. Stability of OSSE spectra was checked using 
hardness ratios (154-282/52-154 keV) for per-orbit spectra being 
added. If orbital spectra within 12 hours were apparently different, 
we lowered the total integration time to include only similar data. 

We used PCA data from 3 to 20 keV, HEXTE data from 20- 
25 to 200 keV and OSSE data from 50 to 1000 keV. In addition, 
four GingafLAC and OSSE simultaneous observations from 1991, 
previously studied bv lGierlinski et al] 1 19971) . were also analyzed. 
The GingafLAC data are available from 1.7 keV, but we decided to 
use exactly the same energy interval as covered by the PCA. For 
the spectral analysis, we use XSPEC 11.3.1k <ArnauJl99^) . 

The spectra may be separated into two groups based on the 
difference of their photon index (see Fig.0. Spectra from 1991 
and 1997, represented by observation 6, have V ~ 1.6 (usual for 
the hard state) and ones from 1999 (observations 25 and 31) have 
r ~ 1.8-2.2, and hereafter we call the corresponding state flat 
(because the spectrum is nearly flat in the EFe plot). The spectra 
from 1996 and 1998 are close to the hardest ones from 1999. We 
do not consider here the broad-band soft state spe ctra, for which 
the sim ilar analysis has been performed by G99 and F rontera et alJ 
1200 ll) . 
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Table 1. Observation log 



No 


RXTE 


Ginga or PCA 


HEXTE 


Date 


Time [UT] 


CGRO 


OSSE 


Time [UT] 




ObsID 


exposure" [s] 


Exposure 6 [s 


1 




VP 


exposure [s] 












1991 hard state 










Gl 




2304 




June 6 


00:18-02:10 


02 


4121 


00:03-02:11 


G2 




888 




June 6 


04:43-06:29 


02 


4040 


04:29-06:51 


G3 




2828 




June 6 


11:03-14:25 


02 


5975 


10:43-14:32 


G4 




1272 




June 6 


20:22-20:44 


02 


1629 


20:02-20:33 










1996 flat state 










1 


10238-01-08-000 


14592 


5083 


March 26 


10:12-17:36 


516.5 


17538 


07:22-20:01 


2 


10238-01-07-000 


8446 


2922 


March 27 


23:06-05:20 


516.5 


23018 


20:20-07:47 


3 


10238-01-06-00 


11455 


3713 


March 29 


11:43-17:33 


516.5 


21868 


08:24-21:06 


4 


10238-01-05-000 


10397 


3524 


March 30 


19:54-01:58 c 


516.5 


24153 


16:51-04:12 c 










1997 hard state 










5 


10239-01-01-00 


9095 


- 


Feb 2 


20:13-02:03 


612.5 


47864 


16:53-04:52 c 


6 


10238-01-03-00 


6441 


1938 


Feb 3 


19:30-22:06 


612.5 


44530 


14:38-02:36 c 


7 


30158-01-01-00 


1175 


809 


Dec 10 


07:08-08:30 


705 


7552 


02:55-12:46 


8 


30158-01-02-00 


2012 


823 


Dec 11 


07:06-08:45 


705 


12004 


04:11-14:05 


9 


30158-01-03-00 


2027 


706 


Dec 14 


08:48-10:20 


705 


11599 


03:09-11:35 


10 


30158-01-05-00 


2614 


901 


Dec 15 


05:26-07:09 


705 


16721 


23:3 1—1 1:18° 


11 


30158-01-06-00 


3210 


941 


Dec 17 


00:40-02:05 


706 


13017 


20:13-07:32 d 


12 


30157-01-02-00 


2309 


784 


Dec 18 


07:07-08:16 


706 


9924 


01:49-08:46 


13 


30158-01-07-00 


2275 


766 


Dec 20 


07:11-08:29 


706 


11847 


02:46-12:52 


14 


30158-01-08-00 


2581 


878 


Dec 21 


05:28-07:05 


706 


14924 


00:51-12:31° 


15 


30157-01-03-00 


2846 


860 


Dec 24 


21:24-23:03 


707 


10127 


21:36-03:18 c 


16 


30161-01-01-000 


13244 


4136 


Dec 28 


13:56-21:03 


707 


13717 


ll:33-00:24 c 


17 


30158-01-12-00 


2836 


916 


Dec 30 


03:52-05:00 


707 


8979 


23:49-07:57 d 










1998 flat state 










18 


30155-01-01-020 


10027 


3897 


Dec 23 


00:07-05:58 


804 


2957 


21:52-09:00 d 


19 


30155-01-02-00 


9625 


3399 


Dec 28 


01:40-07:09 


804 


5124 


23:14-10:29 d 


20 


30161-01-03-01 


9863 


3371 


Dec 28 


13:08-18:37 


804 


5289 


11:50-21:31 


21 


40100-01-04-00 


8472 




Dec 29 


01:37-05:57 


804 


5476 


22:51-10:()8 d 


22 


40100-01-05-00 


8476 




Dec 30 


01:38-05:57 


804 


5288 


22:30-09:46 d 


23 


40100-01-06-00 


7245 




Dec 31 


03:12-07:11 


804 


5552 


23:44-1 l:00 d 










1999 flat state 










24 


40101-01-09-00 


2405 


665 


Oct 5 


18:39-19:45 


831.5 


29605 


15:04-01:12° 


25 


40101-01-11-00 


731 


165 


Oct 6 


19:22-20:03 


831.5 


25960 


14:39-00:50 c 


26 


40101-01-12-00 


865 


- 


Oct 7 


07:10-07:40 


831.5 


28389 


03:08-12:55 


27 


40101-01-15-00 


741 


218 


Oct 8 


08:07-08:48 


831.5 


31014 


04:20-14:35 


28 


4U1U1-U1-10-UU 


757 




Oct 9 


HO. A 1 1 A-Tl 

Uy.41— lU.Zz 


831. j 


33 / l\) 


03.33—13.48 


29 


40099-01-20-01 


1228 


354 


Oct 12 


17:33-19:56 


831.5 


37186 


12:ll-01:23 c 


30 


40100-01-11-01 


4300 




Oct 28 


10:46-15:12 


832 


19559 


06:16-17:46 


31 


40099-01-22-00 


1444 


520 


Nov 8 


14:44-15:25 


832 


17504 


09:52-21:28 


32 


40099-01-23-01 


3084 


1453 


Nov 23 


15:26-17:32 


834.5 


8734 


15:14-21:54 


33 


40100-01-13-01 


770 




Nov 24 


20:08-20:59 


834.5 


19652 


14:51-02:13° 


34 


40100-01-14-02 


479 




Nov 25 


20:05-20:56 


834.5 


16683 


14:24-01:47° 


35 


40100-01-15-03 


1729 




Nov 26 


21:40-22:32 


834.5 


24372 


17:11-04:33° 


36 


40100-01-16-02 


1448 




Nov 27 


19:59-20:50 


834.5 


19063 


15:13-02:34° 


37 


40100-01-17-03 


1902 




Nov 28 


21:32-22:27 


834.5 


10309 


19:19-00:28° 


38 


40100-01-18-03 


1961 




Nov 29 


21:29-22:25 


834.5 


19999 


16:03-03:19° 



a The deadtime-corrected exposure. 

b HEXTE cluster exposure. Symbol '-' means that no HEXTE data are available. 
° The observation finished on the following day. 
d The observation started on the previous day. 
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Figure 2. Correlations between parameters for the model (pexrav). (a) The reflection amplitude, R, vs. the photon spectral index, T. (b) The relativistic 
smearing Gaussian width cr vs. reflection R. (c) The Fe Ka 6.4 keV line equivalent width EW vs. R. 



3 SPECTRAL ANALYSIS 

In order to describe the observed broad-band spectra we use the 
following models: 

(0) Power-law (without or w ith an ex ponential cuto ff) an d 
Compton reflection (pexrav model, Maadziarz & Zdziarski 1995); 

(1) two thermal Comptonization models and reflection, with 
the soft excess also modelled by thermal Comptonization; 

(2) thermal Comptonization with reflection plus a nonthermal 
Comptonization correspondi ng to the soft state spectrum. 

We use the eqpair code lCoppill9 99: G99) for modeling ther- 
mal and nonthermal Comptonization. All models include also a 
Gaussian line at 6.4 keV and interstellar absorption with column 
density TVh which we find often larger than the value of 0.6 x 10 22 
cm" 2 (derived from the redden ing towards the companion star, see 
iBahiciriska-Church et alll995t) . To avoid unreasonably low values 
of this parameter, its low limit was set to 0.5 x 10 22 cm -2 . The pre- 
sented uncertainties are given at a 90 per cent confidence level for 
a single parameter (A^ 2 = 2.71). Fluxes, unless stated otherwise, 
correspond to the range covering all the model emission. 

3.1 The power-law and reflection model in the 3-20 keV data 

The spectra of Cyg X-l clearly show correlations between reflec- 
tion amplitude an d the spectral index ( GCR99, ZLS99). Detailed 
analysis confirms I Zdziars ki et all2003f) that the extend of correla- 
tion is much larger than typical errors in the best-fitting parameters. 
However, since the PCA spectrum falls in a quite narrow energy 
interval and it is difficult to distinguish between various spectral 
components that may form an "effective" power-law. Therefore, it 
is not certain that the values for T and R obtained from the simple 
power-law/reflection fits indeed correspond to the actual physical 
situation. 

Still, in order to compare our results with those of previ- 
ous analysis, we have performed fits similar to those presented in 
ZLS99 and GCR99. We later compare them to the results obtained 
with more physical models (see Sect. 14.41 . We use the XSPEC 
model phabs*(pexrav+gaussian) (model 0), i.e., a power law 
with the photon index, F, and Compton reflection with the rela- 
tive strength, 7?. I Maadziarz & Zdziarski 1995), accompanied by a 
Gaussian fluorescence Fe Ka line (characterized by the relativistic 
smearing width a and the equivalent width EW), all absorbed by 
interstellar material of column density TVh- Hereafter, we assume 
the disc inclination of i = 50° and neutral reflector. 



The fit results are presented in Table |2| Only the low-energy 
3-20 keV (RXTE/PCA and Ginga/L AC) data were fitted, and there- 
fore we did not apply any high-energy cutoff to the power law. We 
find that the correlations between spectral parameters (see Fig.|2j 
are similar to GCR99 results. 

As shown in Sect. 13.21 below, there is a likely overlap of dif- 
ferent spectral components in the PCA energy range. Therefore, 
the model cannot represent a good approach for physical inter- 
pretations for the spectra. Taking into account also the OSSE data 
and including an exponential cutoff in the model, we can model the 
joint data only very roughly, with x 2 /dof ~ 2. This is likely to be 
due to the the shape of the exponential cutoff (which is assumed 
in pexrav) being substantially differen t from the shape of th e cut- 
off of thermal Comptonization (see e.g. lZdziarski et all2003l) . This 
provides an argument against utilizing simple phenomenological 
models in the analysis of BHB spectra. 

3.2 Comptonization model and broad-band spectra 

The hard state spectra are well described by thermal Comptoniza- 
tion iGierlihski et"aT]ll997t FoutanerJll99a iFrontera et alJl200ll) . 
with a weak soft excess. To describe Com ptonization we use the 
XSPEC model eqpair (see lCoppilfl999T: G99). The spectrum of 
seed photons is from a pseudo-Newtonian accretion disc, see G99. 
Parameters of emission are expressed through the compactness, 



where L is the source luminosity, 1Z is the radius of the emitting 
spherical cloud, and err is the Thomson cross section. We consider 
here thermal plasma. The model is characterized by the following 
parameters: t s - the compactness of soft seed photons (assumed 
here to be 1); £h/is - the ratio of the dimensionless energy dissi- 
pation rate in a hot cloud to £ s ; fcT max - the maximal colour tem- 
perature of the disc (at 9.5 gravitational radii, R g — GM/c 2 , or 
at the inner radius of the disc, if it is cut off at larger value); r p - 
the Thomson optical depth corresponding to the ions; R - the re- 
flection amplitude. Since there are large (~25% at la) errors on 
fcT max , we decided to fix it at the common value of 200 eV. The 
model computes the coronal temperature, kT c , and the total opti- 
cal depth, r, from the background electrons and produced e pairs 
self-consistently from the energy and pair balance. For the assumed 
compactness, pair production is negligible for all considered spec- 
tra, and thus the resulting total optical depth r = r p . In the fits, we 
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Table 2. The best-fit parameters for model (fitted to the low-energy Ginga and RXTE/PCA data only). 



No. 


TVh [10 22 cm" 2 ] 


r 


R 


a [keV] 


EW [eV] 


X 2 /dof 


Gl 

G2 
G3 
G4 


0.5t ' 2 
0.5t ' 2 

1-315:1 

1 4+ 12 


1 co+0.09 
l.UO_ 03 

-i R9 +0.09 
1 - oz -0.03 

1 60 +0 10 

J - DU -0.07 
1 RQ+0.09 


1991 

S9+ 019 

u - oy -0.08 

S7+ 018 

uo ' -0.09 

26+ 17 
41 + 019 
1996 


0.46«;« 
0.58«;i 
60+ ' 95 
42+ 97 


152l 5 7 | 
1421 77 

on+36 


3/11 
4/11 
3/11 
5/12 


01 


1.2 ± 0.2 


1.80 ± 0.02 


0.51 ±0.05 


0.68 ±0.16 


96 ± 19 


24/43 


02 


0.6 ±0.2 


1.79 ±0.02 


0.48 ± 0.05 


0.80 ±0.12 


148 ± 20 


20/43 


03 


0.7 ±0.2 


1.80 ±0.02 


0.48 ± 0.04 


0.72 ± 0.12 


138 ± 19 


20/43 


04 


1.2 ±0.2 


1.89 ±0.02 


0.66 ±0.05 
1997 


0.87 ±0.16 


112 ± 19 


24/43 


05 


1.3 ±0.3 


1.65 ±0.02 


0.25 ±0.04 


o 4n +0 - 22 

u -^ u -0.40 


66 ±20 


19/39 


06 


1.2 ±0.3 


1.65 ±0.02 


0.24 ±0.04 


0-381°;^ 


76 ±20 


27/39 


07 


1.0 ±0.3 


1.71 ±0.02 


0.35 ±0.05 


^o +0 - 21 

u - dz -0.26 


79 ±22 


16/39 


08 


1.0 ±0.3 


1.69 ±0.02 


0.29 ±0.04 


n 44+ 019 

u - 4y: -0.24 


86 ±21 


35/39 


09 


1.2 ±0.3 


1.71 ±0.02 


0.36 ±0.05 


56+ 21 

U.OO_ 25 


82 ±22 


26/39 


10 


1.2 ±0.2 


1.69 ±0.02 


0.32 ±0.05 


43+0 20 

u -^°-0.27 


77 ±20 


22/39 


11 


1.2 ±0.2 


1.70 ±0.02 


0.32 ±0.04 


S5+0- 20 


69 ± 19 


15/39 


12 


3.1 ±0.2 


1.72 ±0.02 


0.38 ±0.05 


n+0.42 
U -0 


32 ± 12 


23/39 


13 


2.1 ±0.3 


1.71 ±0.02 


0.37 ±0.05 


u - ou -0.30 


56 ± 19 


21/39 


14 


1.1 ±0.3 


1.70 ±0.02 


0.37 ±0.05 


43+0-20 

u -^°-0.29 


78 ±21 


21/39 


15 


2.1 ±0.3 


1.70 ±0.02 


0.35 ±0.05 


i^+°- 53 
u.xo_ 15 


41 ± 18 


17/39 


16 


2.1 ±0.2 


1.69 ±0.02 


0.41 ±0.05 


17+0 43 


57 ± 16 


22/39 


17 


2.7 ±0.3 


1.71 ±0.02 


0.36 ±0.05 
1998 


03+ 2 79 

u - uo -0.03 


37 ± 17 


19/39 


18 


1.2 ±0.2 


1.87 ±0.02 


0.54 ±0.05 


0.73 ±0.16 


99 ±20 


19/39 


19 


1.2 ±0.2 


1.84 ±0.02 


0.48 ±0.05 


0.71 ±0.17 


95 ±20 


19/39 


20 


1.2 ±0.2 


1.85 ±0.02 


0.50 ±0.05 


0.71 ±0.17 


96 ±20 


17/39 


21 


1.1 ±0.2 


1.83 ±0.02 


0.47 ±0.05 


0.71 ±0.17 


93 ±20 


14/39 


22 


0.7 ±0.2 


1.83 ±0.02 


0.47 ±0.05 


0.69 ±0.14 


123 ±21 


20/39 


23 


0.7 ±0.2 


1.82 ±0.02 


0.46 ±0.05 
1999 


0.72 ±0.13 


135 ± 21 


17/39 


24 


0.5+r 


2.09 ±0.01 


0.89 ±0.07 


0.93 ±0.12 


203 ± 23 


28/33 


25 


o.str 


2.30 ±0.02 


1.38 ±0.16 


1.12 ± 0.16 


249 ± 34 


48/33 


26 


o.str 


2.23 ±0.01 


1.15 ±0.10 


1.09 ±0.13 


247 ± 27 


43/33 


27 


0.5+°r 2 


2.05 ±0.01 


0.81 ±0.08 


0.95 ±0.14 


188 ± 25 


14/33 


lb 


n t i n *7 
U. i ± U.o 


A.UO ± U.Uo 


U.84 ± U.1U 


U.oz ± U.lO 


lOl ± ZD 


hi a 


29 


0.51S- 1 


2.01 ±0.01 


0.74 ±0.07 


0.91 ±0.15 


185 ± 25 


22/33 


30 


0.5+T 


2.20 ±0.01 


1.07 ±0.07 


1.03 ± 0.11 


237 ± 22 


27/33 


31 


0.5+°r 2 


1.94 ±0.02 


0.61 ±0.06 


0.82 ±0.13 


170 ± 22 


25/33 


32 


1.1 ±0.3 


1.85 ±0.03 


0.50 ±0.06 


0.54 ±0.23 


79 ±22 


20/33 


33 


0.7t°;| 


1.87 ±0.03 


0.53 ±0.07 


0.63 ±0.19 


113 ± 26 


30/33 


34 


0.61°;? 


1.86 ±0.03 


0.47 ±0.08 


0.61 ±0.19 


125 ± 28 


14/33 


35 


0.5+ , 2 


1.88 ±0.02 


0.52 ±0.05 


0.65 ±0.13 


141 ± 18 


25/33 


36 


0.5^' 2 


1.88 ±0.02 


0.55 ±0.06 


0.72 ± 0.14 


156 ± 22 


25/33 


37 


1.2 ±0.3 


1.92 ±0.03 


0.62 ±0.07 


0.68 ±0.18 


109 ± 24 


27/33 


38 


1.2 ±0.3 


1.98 ±0.03 


0.77 ±0.08 


0.78 ±0.16 


125 ± 24 


39/33 



assume the inner disc radius of 6R g (this parameter has only effect 
on the relativistic smearing of the reflected component, that could 
be hardly resolved with present energy resolution). 

The normalization of eqpair, corresponding to the disc com- 
ponent, is f c M 2 cos i/(D 2 (3c )> where f c is the covering factor, M 
is the black hole mass in units of the solar mass, D is the distance 



to the source in units of kpc, and J3 C is the ratio of the colour tem- 
perature to the effective one. 

At low energies the effect of the interstellar absorption is 
clearly visible. The addition of a additional soft component leads to 
further improvement of t he fits. Such soft exce sses were observed 
in both the hard state by iDi Salvo et alj fcOOll) and in a flat state 
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(similar to that analyzed here) bv lFrontera et alJ 1200 lh in Beppo- 
SAX data. This component, that nature we address below, however, 
is relatively weak in case of 1991 and 1997 (observations G1-G4, 
05-17) hard state data. The spectra observed in 1996, 1998 and 
1999 (observations 1-4, 18-38) appear similar in the overall shape 
to those of hard state, but are significantly softer (see Fig. |2j- The 
soft excess in the flat state is stronger. 

We stress that the requirement of an additional soft excess is 
implied only by the joint PCA/HEXTE/OSSE data, since the PCA 
data cover a too narrow energy range. Even if the actual spectrum 
in the PCA 3-20 keV band is not a power-law but is e.g. concave, 
a good fit with a power law plus reflection (model 0) can be still 
achieved. However, the real strength of Compton reflection can be 
significantly different. 

Since the data require an additional component only in a rela- 
tively narrow, ~3-10 keV, range, the parameters of the soft excess 
cannot be constrained accurately. Below we consider a number of 
physically realistic scenarios of its nature. In each of the considered 
models, we restrict the parameters controlling its spectral shape to 
values that makes its flux significant only at low energies, and we 
fit only its normalization. These fits allow us to completely describe 
the broad-band spectra and to constrain the parameters of main con- 
tinuum. 

3.2.1 High temperature of the optically-thick disc 

The additional component may, in principle, be emitted by the 
hottest part of the optically-thick disc provided its temperature is 
high enough. We find that the spectra of the observ ations 1-4, 18- 
38 can be well fitted with feT max ~ 1 keV. But P i Salvo et"ai] 
1200 lh showed that the spectral decomposition of the BeppoSAX 
data of Cyg X-l in the hard state requires the presence of both 
the soft blackbody disc photons (with fcT max ~ 0.1-0.2 keV, and 
an additional soft excess component. Therefore, we consider this 
model to be not appropriate. 

3.2.2 Hybrid model 

We also tried to apply hybrid fhermal+nonthermal eqpair model to 
the spectra. Most of the data require nearly nonthermal injection 
and the resulting spectrum becomes a power-law directed by PCA 
part of the spectrum. Because the observed hard tail is not power- 
law, large reflection appears to mimic the cutoff region, but it is still 
not enough to describe both the soft excess below 10 keV and the 
hard tail, in which the systematic difference between the data and 
the model remains. Therefore, we ruled out this model. 

3.2.3 Two thermal Comptonization components 

The soft component can be described by additional thermal Comp- 

^^^^^^^^^^^^^^^^^^^ ^ ^^^^^^^n^^^^T, 

tonization I Frontera et al. 2001; Di Salvo et al. 2001). We use the 
model 1, phabs(comptt+eqpair+ga ussian), where e qpair gives 
the main Comptonization and comptt (Titarchuk 1994) - the addi- 
tional soft component. Since the parameters of the soft excess are 
rather weakly constrained by our data, we fixed the parameters of 
comptt at kT c — 20 keV and r = 1. Such a model produces a soft 
power law that does not extend to very high energies. This model 
fitted to the data of the observation 24 is shown in Fig. [5] In spite 
of its simplicity, it yields a relatively good description of the data. 

The fit parameters are given in Table|5| and the dependencies 
between various model parameters are shown in Fig.|4| We also 
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Figure 3. The model spectrum with two thermal Comptonization compo- 
nents (model 1, Sect. 13.2.31 fitted to the RXTE+OSSE observation 24 (from 
1999). The spectral components of the fit are shown by the dotted (green), 
dot-dashed (red), and dashed (blue) curves, which correspond to the main 
thermal-Comptonization continuum, the additional thermal Comptoniza- 
tion with kT c = 20 keV and r = 1, and the Compton reflection including 
the Fe line, respectively. The solid (black) curve shows the total spectrum. 
The lower panel shows the residuals of the fit. 



quote the values of the spectral index V of the power law obtained 
from the least-square fitting of the logarithm of the intrinsic model 
flux at a logarithmic energy grid in the 2-10 keV range (chosen to 
enable comparison with results of other papers). 

The meaning of the normalization of eqpair is described in 
Sect. l3~2l above. Substituting f c = 1, M = 10M^, i = 50°, 
D = 2.0 kpc (see references in G 99, lFrontera et al 1200 lh and f3 c = 
1.7 IShimura & Takaharal ll995). we expect the normalization of 
~1.92. The lower normalization of the obtained fits (see Table[3j, is 
caused either by a larger j3 c ~ 1.77-1.98 or by a smaller covering 
factor f c . A slightly smaller, than the assumed disc temperature 
of 200 eV, can also reduce the normalization. On the other hand, 
some flat state normalizations are larger than the expected value of 
1.92. This cannot be explained by changing f c (smaller covering 
factor can only reduce the normalization), but could be a result of 
somewhat larger inner disc radius or larger temperature. The largest 
observed normalizations correspond to fcT max ~ 220 — 245 eV. 

We see that there is a clear anti-correlation between the 
strength of the soft component and hardness of spectrum expressed 
in terms of £h/£ s . There is also a correlation between r and l-^/ls- 
The reflection amplitude is correlated with the Fe line equivalent 
width EW and anti-correlated with i\y/l s (while at high values of 
l-h/ls the anti-correlation possibly breaks down). 

The electron temperature kT a is a calculated parameter and no 
errors on it can be obtained from fitting. However, we estimated its 
lcr limits using its extremal values within the uncertainties of the 
parameters controlling spectral shape, i.e., £h./l& and r. This esti- 
mation gives us a possible conservative error on kT c of about 15 
keV for both models 1 and 2. Taking this into account, no correla- 
tions between fcT c -flux and kT c — £h jl% are apparent. 
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Table 3. The best-fit parameters for model 1. 
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a Hydrogen column density, in units 10 22 cm~ 2 . 

b Normalization of the eqpair model component. 

c The unabsorbed total model flux, in units of 10 _8 erg cm~ 2 s _1 . 

d The unabsorbed model flux of the comptt component, in units of 10~ 8 erg cm~ 2 s _1 . 

e Temperature of the emitting plasma in keV (for the eqpair component). 

f Photon spectral index of the eqpair component in the 2-10 keV range. 
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Figure 4. Dependencies of the spectral parameters for model 1 (Sect. 13.2. j| fitted to the 1991 and 1997 (blue filled circles) and the 1996, 1998 and 1999 data 
(red filled squares) on Compton amplification factor of the main Comptonization component i^/is. (a) The reflection fraction R; (b) the relativistic smearing 
Gaussian width a at 6.4 keV; (c) ratio of the additional thermal Comptonization flux to the total flux; (d) Thomson optical depth of the main Comptonization 
continuum component r; (e) electron temperature of the main Comptonization component fcT c ; (f) total luminosity of the Comptonizing cloud (assuming 
D =2 kpc). 
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Figure 5. The model spectrum with the thermal and nonthermal Comp- 
tonization components (model 2, Sect. 13.2.41 fitted to the RXTE+OSSE ob- 
servation 24 (from 1999) together with the COMPTEL IMcConnell et alj 
l2002t) hard-state data (marked by squares). The spectral components of the 
fit are shown by the dotted (green), dot-dashed (red), and dashed (blue) 
curves, which correspond to the main thermal-Comptonization continuum, 
the additional nonthermal Comptonization, and the Compton reflection in- 
cluding the Fe line, respectively. The solid curve shows the total spectrum. 
The lower panel shows the residuals of the fit. 



3.2.4 Thermal continuum and nonthermal Comptonization 
component 

As an alternative to the second thermal-Compton component, we 
consider here addition of a nonthermal Comptonization compo- 
nent. We find that such a nonthermal component can describe both 
the soft excess and th e MeV nonthermal tail observed in hard 
states by COMPTEL IMcConnell et alJl2002l) . while neither of 
these components can be described by the main thermal-Compton 
emission. This model fitted to the spectrum 24 is shown in Fig. [5] 
Note that the COMPTEL data are shown for illustration only and 
were not taken into account in the fitting. The used model 2 consists 
of phabs(eqpair+eqpair+gaussian), in which the second eqpair 
component produces the nonthermal spectrum. In eqpair, the avail- 
able power is supplied in part into heating electrons and in part into 
their acceleration, with the resulting steady-state electron distribu- 
tion calculated self-consistently. The compactness corresponding 
to the acceleration is hereafter denoted as £ nt h- Then the relative 
fraction of the input power going into the nonthermal acceleration 
is 4ith/^h, where £h (as before) corresponds to the total rate of 
energy dissipation in the plasma. 

For that additional component, we assumed that all the avail- 
able power goes into nonthermal acceleration, i.e., £ n th/(-h = L 
Note that the resulting self-consistent electron distribution is not 
purely nonthermal but hybrid, i.e., it does contain a low-energy 
Maxwellian heated by Compton and Coulomb interactions. We 
further assume R = 0, the power-law index of the acceler- 
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Figure 6. Dependencies of the spectral parameters for model 2 (Sect. 13.2.41 fitted to the 1991 and 1997 (blue filled circles) and the 1996, 1998 and 1999 data 
(red filled squares) on Compton amplification factor of the main Comptonization component 1^1 is- The meaning of the axes is the same as in Fig.|4] except 
(c) ratio of the additional nonthermal Comptonization flux in model 2 to the total flux. 



ated electrons i of = 2.4 (see iMcConnell etai]l2002t G99; 
iFrontera et aljEoOll IPoutanen & Coppilll998l) . the minimum and 
maximum Lorentz factors of the power law of / ymin — 1.3 and 
7max = 1000, respectively, fcT max equal to that of the main com- 
ponent, £h/£s = 1, and r p = 1. No pair production is required, and 
t is found to be equal to r p . 

The best-fit parameters are presented in Tableland correla- 
tions between them are shown in Figs.lrjland ll II The normalization 
of the thermal eqpair component corresponds to f3 c = 1.63-2.01, 
the values similar to those obtained with model 1. The strength 
of the additional (non-thermal) component is again clearly anti- 
correlated with the spectral hardness. The anti-correlation between 
R and lh/l s is apparent and corresponds to the R — V correlation 
that we discuss in Sect. 14.41 The values of R and the Fe line equiv- 
alent width EW are well correlated and can be approximated by a 
linear function passing through zero (see Fig. ll lb ). Again, there is a 
correlation between r and £h/^s- The electron temperature does not 
seemingly depend on the hardness, but the spread becomes smaller 
at larger 4/4 where kT c ~ 90 keV. 



4 DISCUSSION 
4.1 Absorption 

As noted in Sect. [3] the hydrogen column density JVh was free in 
our fits. We find that our data require absorption significantly larger 
than 0.6 ± 0.2 x 10 22 cm" 2 which is derived from the red dening 
towards the companion star (Balucinsk a-Church et alll 995). 




Figure 7. (a) Fitted value of absorption column density iVfj for model 1 
vs. orbital phase of the binary system, (b) Same, for model 2. Distributions 
are repeated twice for clarity. iVjj is in units of fO 22 cm -2 . 
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Table 4. The best-fit parameters for model 2. 
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a Hydrogen column density, in units 10 22 cm -2 . 

b Normalization of the thermal eqpair model component. 

c The unabsorbed total model flux (without reflection) in units of 10~ 8 erg cm -2 s — 

d The unabsorbed model flux from the nonthermal component, in units of 10 — 8 erg cm -2 s — 1 . 

e Temperature of the emitting plasma in keV (for the thermal eqpair component). 

' Photon spectral index of the thermal eqpair component in the 2-10 keV range. 



Cyg X-l is known to show X-ray dips in its light curve, caused 
by obscuration by the stellar wind from the companion star. Dur- 
ing dips the absorption increases up to A^h = 20 x 10 22 cm -2 . 
The absorption shows strong orbital phase dependence with max- 
imal column density around phase 0, when companion star is in 



front of the black hole <Feng & Cul200llBaiucinska-Church et alJ 
2000 ). We plotted the fitted values of TVh versus Cyg X-l orbital 
phase (see Fig. Q, using the ephemeris with the time of the pri- 
mary minimum at 50234.79 MJD and the p eriod of 5.599829 days 
jBrocksoPD et alJll999tlLaSala et alj|l 9-98). We see a good corre- 
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lation of TVh with the orbital phase, which indicates that variable 
absorption indeed can be caused by the companion's wind obscur- 
ing the X-rays from the black hole. 

We note that since our models are relatively complicated in 
the range of 3-10 keV (consisting of absorption, Comptonization 
together with the seed photon emission and a soft excess), and also, 
the PCA energy range is affected by absorption rather weakly, it 
is difficult to determine the exact values of Nh- Even though, the 
derived values of the parameter are in the range quoted by other 
authors and its relative changes are quite remarkable. 



4.2 Spectral variability patterns 

Using CGROfBKTSE and RXTEl ASM data. lZdziarski et aljj2002[) 
showed that Cyg X-l has two types of variability - changes of flux 
without spectral slope change and pivoting at ~ 50 keV which pro- 
duces anti-correlation of the fluxes in the soft and hard part of the 
spectrum. However, those instruments do not provide detailed spec- 
tral information giving fluxes only in some energy intervals. The 
effective photon spectral indices using fluxes in two energy bands 
and calculated in the wide (20-300 keV) energy interval, where real 
spectra experience a cutoff, should be tre ated with caution. Now we 
have a possibility to check the results of Zdziarski et al. 1 200^) us- 
ing our set of observations, on which we have a detailed spectral 
information in the wide energy range. For this purpose we are us- 
ing the model spectra obtained from model 1 . 

In Fig. [3] we present the spectra related to different time peri- 
ods. It is possible to see that 1991, 1996 - 1998 spectra only change 
their normalization. However, spectral slope is different for various 
years and forms two groups: 1991+1997 spectra and 1996+1998 
ones. In the 1999 data, we again can see normalization changes 
(see pairs of red solid curves) as well as pivoting behaviour. 

The T-flux correlations are shown on Fig. [5] It is clearly seen 
that the 1991 and 1996-1998 data do not show dependence between 
flux and F. The 1999 data show clear anti-correlation between F 
and flux on low energies (3-12 keV) and correlation on high ener- 
gies (20-100, 100-300 keV), that indicates the pivoting behaviour 
with the piv ot energy between 12 and ~ 50 keV. On Fig. |5Jl (see 
also fig. 8 in lZdziarski et al.l2002l) . it can be seen that spectra from 
1996 and 1998 are somewhat softer at high energies than the hard 
state spectra from 1991 and 1997, but the fluxes in all energy in- 
tervals generally correlate with each other. The 1999 spectra show 
clearly different dependence: the high energy fluxes are correlated, 
but there is a clear anti-correlation between 3-12 keV and 100-300 
keV fluxes due to pivoting. The timing behaviour has also changed 
its nat ure in late 1998, as was pointed out by Pottschmi dt et alJ 
120031) . 

4.3 QPO frequencies 

Among the data sets we have studied, there are timing 
data for 33 observations avail a ble fr om the recent paper by 
lAxelsson. Borgonovo & Larssonl l2005t) . This allows us to check 
the relation between the characteristic frequencies in the power 
densit y spectra and the spectral parameters. In the work of 
lAxelsson et alJ {2005), several values of QPO frequencies ^qpo 
might be determined during one observation, and for these data 
points we assume the averaged middle frequency and consider the 
uncertainty from lower to higher of obtained values. In agreement 
with earlier results of GCR99, we find a clear anti-correlation be- 
tween the characteristic frequencies of the aperiodic variability and 
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Figure 8. Sample of the spectra from different periods. Lowest and highest 
spectra from each time period are shown. Red solid lines - 1999; blue dotted 
lines - 1991 (the 1997 spectra are similar to them); black dashed lines - 
1996 and 1998. The 1991+1997 and 1996+1998 data form two groups with 
slightly different spectral slopes. 

l\i/is (Fig. I10> . indicating an apparent connection between QPO 
frequencies and the parameters of the Comptonizing region. This 
provides an argument (but not a proof) in favour of the presence of 
a hot inner corona and a variable inner radius of the surrounding 
disc. 

The frequency-hardness correlation can be described by a 
power-law vqpo oc (4/4)~ a with a = -1.48 ± 0.04. The best 
fit is shown by the solid curve on Fig . Hol Simil ar correlation was 
observed bv lPottschmidt et alji2003l) and CMowak et alJl2002h . 

4.4 Comparison between phenomenological and physical 
spectral models and R — F correlation 

Fig. II ll compares the results obtained with the simple phenomeno- 
logical model of power law + reflection in the 3-20 keV energy 
range (model 0) with those from our model 2 applied to the 3- 
1000 keV range. More elaborate and physically justified models 
utilizing the full energy range of our data do not change the pic- 
ture qualitatively. On the quantitative level we find that the simple 
power law + reflection spectral fits to the 3-20 keV data overesti- 
mated the amplitude of the reflected component R and the slope T 
of the primary Comptonization continuum. We confirm, however, 
that the simple models did rank correctly the spectra according to 
the strength of the reflected component and slope of the Comp- 
tonized radiation, as it was demonstrated in the original publica- 
tions on this subject (ZLS99; GCR99). It is also illustrated by the 
lower three panels of Fig. 1111 The difference of the obtained pa- 
rameters comes from the fact, that for wide-energy observations 
the main thermal Comptonization component that describes well 
the hard energy tail may lie well below the observed flux in the 2- 
10 keV range (see Fig. [3] and |5J and has a different slope in this 
band. The difference is largest for 1996, 1998 and 1999 data, while 
in the 1991 and 1997 cases the soft excess is weak and parame- 
ters obtained with the physical models are similar to those obtained 
with the phenomenological ones. 

On Fig. [TTJt, we see a clear correlation between F and R. 
For comparison, we also show the dependencies predicted by 
the plasma ejection model of Beloborodov 1 1999a b) and MBP01 
(cylindrical geometry with h/r = 2), using the geometric param- 
eter of that model of fi s — 0.4 and 0.5, i = 50°, albedo of the 
reflecting medium of 0.15 and r = 2. We used the dependence 
between the amplification factor A of Comptonization and F from 
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Figure 9. Flux-spectral index (3-12 keV) and flux-flux correlations. Respective energy intervals are indicated along the axis. The 1991 and 1997 data indicated 
by blue rilled circles, 1996 and 1998 - by red filled squares and 1999 - by red open squares. All fluxes are in units of keV cm -2 s — 1 , 
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Figure 10. Dependence between QPO frequencies and £h/l s < The i^/is 
are best fit v a lues f rom model 2. The QPO frequencies are from 
lAxelsson et alJ l2005t) . The solid curve shows the best power-law fit 

^qpo oc (A/4)- 1 - 48 . 



MBP01. We also compared our data with the dependence expected 
in the model of ZLS99, assuming the black body temperature of 
0.2 keV, appropriate for Cyg X-l, with one minor change. In the 
original paper, all the reflection luminosity was assumed to reach 
the observer. Reflection amplitude is an integral that consists of 
two parts, from the disc inside the corona and from the outer part 
of it. We multiplied the part of luminosity coming from the former 
part by e~ T , to approximately take into account scattering of radi- 
ation in the corona (r = 1 was chosen). We see that this model 
cannot quantitatively describe the presented data. Moreover, taking 
into account intrinsic dissipation in the disc (see Appendix IaI for 
details) will further increase the slope of the dependence making 
the discrepancy larger. Intrinsic dissipation becomes important for 
a small inner disc radius (when reflection is relatively large) and the 
increase of soft seed photon flux in that case makes the spectrum 
softer (see Beloborodov 2001). 



4.5 Physical scenario 

The hard spectral state of black hole binaries is commonly de- 
fined as the state in which the spectrum is dominated by the 
hard Comptonization component, without significant contribution 
of the blackbody-type emission from the optically thick accre- 
tion disc. Naturally, the hard state is not characterized by a single, 
uniquely defined spectrum, but rather includes a continuum of spec- 
tral shapes with the major spectral parameters vary in a rather broad 
range. The diversity of the hard state spectra reflects the position of 
the source with respect to the "bottom hard" state and the soft state. 
Quantitatively, this position can be characterized by the strength 



14 A. Ibragimov et al. 




T (model 0) R (model 0) EW (eV, model 0) 

Figure 11. (a-c) Correlations obtained using model (green open circles), model 1 (blue filled circles) and model 2 (red filled squares), (a) The spectral slope 
r vs. reflection scaling factor R correlation. For model 0, T is a fitting parameter, for model 2 - index of a power-law fitted to the spectral model in the 2-10 
keV range. The solid curve shows the dependence predicted for the ejection model I Beloborodov 1999a; MBP01) with the parameters of i = 50°, the albedo 
of 0.15, fis = 0.5 and t = 2, the dotted curve - the same for fi B = 0.4. The dashed curve shows the dependence predicted by the model with a hot spherical 
corona and cold overlapping disc (Pouta nen et aljl997t ZLS99) with the black body temperature of 0.2 keV, see Sect. I4.4l for details; the dot-dashed curve 
shows the same model with the dissipation effect taken into account (dissipation parameter Li nt = 1, Li nt /47r = 1 at r < f , see AppendixlAlfor details), (b) 
The relativistic smearing Gaussian width a at 6.4 keV vs. R. (c) The equivalent width of the 6.4 keV line EW vs. R. The straight line is Blf [eV] = 500i?. 
(d-f) r, R and EW obtained from model 2 vs. those from model 0. 



of the reflected component (reflection scaling factor R or Fe line 
equivalent width EW) or properties of the main Comptonized 
component (Comptonization parameter, or Compton amplification 
factor A = £h/l$, or the photon index T in the low energy limit) 
or characteristic frequencies of the aperiodic variability. Existence 
of good correlations between all these quantities suggests that they 
all are an equally good measure of the source position within the 
hard state. The results found in earlier work (ZLS99; G CR99; 
iGilfanov et aljEoOol iGilfanov. Churazov & Revnivtsevll2004) and 
presented in the previous sections of this paper suggest the fol- 
lowing pattern of spectral and temporal variability. Increase of the 
strength of the reflected component is accompanied by the increase 
of the width of the Fe line, increase of the characteristic QPO fre- 
quencies and softening of the Comptonized component observed as 
increase of its photon index F in the low energy limit. 

We find from our spectral analysis that in the "bottom hard" 
state the broad band spectrum (3-1000 keV) is mostly described 
by the single thermal Comptonization spectrum with superimposed 
component due to reflection of the primary emission from relatively 
cool and neutral, or partly ionized, optically thick matter (the ac- 
cretion disc), with an additional relatively weak soft component. 
As the source moves towards the soft state, the strength of the 
reflected component increases, and the soft component becomes 
more significant. Considering the 3-1000 keV energy range cov- 



ered by our data, this spectral component reveals itself most clearly 
in the E < 10 keV energy domain as the "soft excess". Another 
independent indication of spectral complexity is the gamma-ray 
power -law tail detected at Me V energies by the COMPTEL tele- 
scope iMcConnell et all2002t) . 

From the point of view of the formal fit quality, the E < 10 
keV excess can be described equally well by an additional thermal 
Comptonization component due to low temperature, low Comp- 
tonization parameter plasma or by non-thermal Comptonization 
with the power law index of accelerated electrons Tinj ~ 2 — 3. Ow- 
ing to the complex shape of the continuum at these energies defined 
by the superposition of several spectral components, these two pos- 
sibilities can not be easily discriminated based solely on the low 
energy data. The task is further complicated by the rather limited 
low energy coverage provided by the PCA instrument, E > 3 keV. 
However the above possibilities predict very different behaviour in 
the ~ MeV energy domain, where the main thermal Comptoniza- 
tion component diminishes and the power law tail due to the non- 
thermal Comptonization should reveal itself. As the OSSE sensitiv- 
ity and energy range are insufficient to probe existence of the MeV 
tail correlated with the E < 10 keV excess, we can not, strictly 
speaking, give preference to either of these two models. 

There are however several additional considerations to be 
taken into account: (1) COMPTEL detected a weak MeV tail in 
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the averaged hard state data for Cyg X-l iMcConnell et all2002l) . 
The slope and amplitude of this tail is qualitatively consistent with 
the extrapolation of the non-thermal Comptonized component, re- 
quired to explain the E < 10 keV excess (Fig[5J; (2) in the soft state 
the non-thermal power law is the dominant (the only) high energy 
component; (3) the E < 10 keV excess is more pronounced in the 
spectra characterized by large reflection and rather steep slope of 
the main Comptonized component, i.e. in those sufficiently close 
to the soft state. Its strength seems to increase with increase of the 
reflection. These arguments suggest that the non-thermal origin of 
the E< 10 keV excess is more plausible. We note, that the 1991 
and 1997 data show much weaker excess. This may be explained 
by the lower relative luminosity of the non-thermal Comptonization 
component which therefore reveals itself at lower energies, below 
the E — 3 keV threshold of the PCA instrument, but can be de- 
tected by instr uments which have response at lower ene rgies, i.e. 
BeppoSAX (see lFrontera et alj200lHDi Salvo et all200ll) . 

The overall qualitative picture can be outlined as follows. 
The overall geometry of the accretion flow is adequately repre- 
sented by the truncated disc model with the inner radius of the 
standard optically thick geometrically thin disc varying from ~ 3 
R g to ~ several tens R g . Inside this radius the accretion flow 
proceeds via quasi-spherical optically thin hot flow. The plausi- 
ble mechanism governing the transition from the disc accretion to 
the coronal flow is the disc evaporation process as proposed by 
iMever & Mever- Hofmeisterl ll994t) . The geometrically thin disc 
gives rise to the soft black body type component. In addition, due to 
dynamo, sol ar-typ e magnetic flares can be produced above the ac- 
cretion disc I Gal eev et all 19791) . The electrons there can be acceler- 
ated and form non-thermal distribution. Comptonization of the disc 
emission on these electrons results in the power law-like Comp- 
tonized emission. The inner optically thin flow gives rise to the 
thermal Comptonization component. The relative contributions of 
non-thermal and thermal Comptonized components are defined by 
the fractions of the gravitational energy released in the disc (i.e. 
outside Ri n ) and in the inner hot flow (inside Ri n ). The position 
of the transition radius is defined by the mass accretion rate and 
is modified by the irradiation-related effects. The transition radius 
decreases as the mass accretion rate increases. 

The QPOs are due to some processes in the transition region 
near R[ n and approximately scale with the Keplerian frequency and 
other characteristic time scales of the coronal flow and standard 
accretion disc in the transition region. 

The configuration with the large inner disc radius, probably 
_Ri n > 50 — 1007? 9 , corresponds to the classical hard state. The 
main features of this "bottom hard" state are low strength of the 
reflected continuum, relatively narrow fluorescent Fe line of small 
equivalent width, large Comptonization parameter of the thermal 
Comptonized component (hard spectra with the low energy photon 
index V ~ 1.6), low frequencies of QPOs. As only small fraction of 
the gravitational energy is released in the disc, the contribution of 
the non-thermal component is small and the spectrum is adequately 
described by thermal Comptonization. 

As the mass accretion rate increases, the transition radius de- 
creases, the disc moves towards the compact object. This results in 
increase of the reflection, broader fluorescent Fe line, larger QPO 
frequencies, smaller £h/£ s , i.e. smaller Comptonization parameter 
in the inner hot flow. The contribution of the non-thermal compo- 
nent increases. The optical depth of the thermal plasma of the inner 
flow decreases due the shrinking of the inner hot flow as the disc 
extends towards the compact object. 

The classical soft state (we ignore all the complications and 



sub-states here) corresponds to the accretion disc extending all the 
way towards the last stable orbit or very close to it. Correspond- 
ingly the inner hot flow disappears and the dominant or the only 
hard component is the one due to non-thermal Comptonization of 
the disc emission on the non-thermal electrons accelerated in the 
magnetic loops/flares above the disc. 

The behaviour of the temperature of the thermal Comptoniza- 
tion component is unclear. It seems relatively constant, which sug- 
gests of possible presence of electron-positron pairs (see MBP01). 

The physical scenario qualitatively outlined above is based on 
the truncated disc picture and on the assumption that the spectral 
evolution is governed by the change of the transition radius be- 
tween the standard accretion disc and the hot inner flow. In this 
picture many of the observed correlations can be explained natu- 
rally. However, the R — V correlation is significantly better quan- 
titatively explained by the non-stationary corona model (MBP01), 
in which the governing parameter is the velocity of the blobs of 
emitting plasma relative to the accretion disc. We note that both 
models are geometrical in their nature and the predicted qualita- 
tive relations between the physical parameters are obtained with a 
number of simplifying assumptions. Therefore results of quantita- 
tive comparison of the model predictions with the observed pattern 
of the spectral variability should be interpreted with caution and 
any conclusions regarding validity of either model based on such a 
comparison would be premature. 



5 CONCLUSIONS 

Based on the broad band (3-1000 keV) data from simultaneous 
observations by Ginga and CGRO/OSSE in 1991 and PCA and 
HEXTE instruments aboard RXTE and OSSE in 1996-1999 we 
studied the spectral variability of Cyg X-l. 

(i) We confirm earlier results on R — T correlation. Consider- 
ing the 3-20 keV data we find very tight one-parameter relations 
between reflection, spectral index and the width of the Fe line. 

(ii) More elaborate and physically justified models utilizing the 
full energy range of our data do not change the picture qualita- 
tively. On the quantitative level we find that the simple power law 
+ reflection spectral fits to the 3-20 keV data overestimated the am- 
plitude of the reflected component 7? and the slope T of the primary 
Comptonization continuum. We confirm, however, that the simple 
models did rank correctly the spectra according to the strength of 
the reflected component and slope of the Comptonized radiation, 
as it was demonstrated in the original publications on this subject 
(ZLS99; GCR99). 

(iii) Based on the analysis of the broad band data we found that 
the spectra in our sample can be adequately described by the ther- 
mal Comptonized component with superposed reflection from the 
optically thick disc and a soft excess. This excess is relatively weak 
in case of hardest spectra of our sample (T ~ 1.7). As the strength 
of the reflection increases, the excess becomes much more signifi- 
cant. Presence of this excess was the primary reason for the simple 
spectral approximations of the 3-20 keV data to overestimate both 
R and T. The nature of this excess cannot be unambiguously de- 
termined from our data. Based on the circumstantial evidence we 
suggest that it is the lower energy part of the non-thermal Comp- 
tonized component with the power law index of accelerated elec- 
trons Tinj ~ 2 — 3. At higher energies this non-thermal component 
reveals itself as a power law detected by COMPTEL at MeV ener- 
gies in the average hard state spectrum of Cyg X-l. 
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(iv) We note the variability of the absorption correlated with the 
phase of binary system. These results confirm previous findings of 
the X-ray dips in the source. 

(v) The overall pattern of spectral and temporal variability can 
be summarized as follows. Increase of the strength of the reflected 
component is accompanied by the increase of the width of the Fe 
line, increase of the characteristic QPO frequencies and softening 
of the Comptonized component observed as increase of its photon 
index V in the low energy limit or, equivalently, decrease of the 
Compton amplification factor t^/la. Simultaneously, the optical 
depth of the thermal Comptonization decreases and the fractional 
contribution of the non-thermal component to the total energy flux 
increases. The exact behaviour of the electron temperature in the 
hot inner flow is not constrained by our data. 

(vi) We suggest a qualitative physical scenario naturally ex- 
plaining the observed behaviour. In this scenario the evolution of 
the spectral parameters is governed by the value of the transition 
radius between the standard optically thick accretion disc and the 
inner quasi-spherical hot flow. The thermal Comptonized compo- 
nent originates in the inner hot flow as a result of Comptonization 
of the soft photons emitted by the accretion disc. The origin of the 
non-thermal component is related to the optically thick disc, for ex- 
ample it can be produced due to non-thermal electrons accelerated 
near the surface of the optically thick disc in the solar-type mag- 
netic flares. The relative contributions of non-thermal and thermal 
components to the total energy flux depends on the fractions of the 
gravitational energy of accreting matter released in the optically 
thick disc and in the hot inner flow. 
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APPENDIX A: DISSIPATION AND ATTENUATION IN 
THE DISC-HOT FLOW MODEL 

ZLS99 have considered an idealized geometrical model for thermal 
Comptonization, reprocessing and reflection in an accretion flow 
consisting of a central hot sphere surrounded by a flat cold disc, see 
fig. 2 in ZLS99. The sphere has a unit radius, and the inner radius 
of the disc can assume any value, d. For d < 1, there is an over- 
lap between the two components. The hot sphere Comptonizes soft 
seed photons emitted by the disc. In the original model of ZLS99, 
the disc reprocesses and reemits only the photons emitted by the 
sphere incident on the disc. 

Here, we generalize that model to include intrinsic dissipation 
in the cold disc (as expected in an accretion flow). Also, we take 
into account scattering of the Compton-refiected photons in the hot 
sphere, which was neglected in ZLS99. For completeness, we give 
here the full set of relevant equations, but refer the reader to ZLS99 
for details of the derivation. 

The hot sphere has a unit luminosity and emits isotropically. 
The total flux incident on the disc at a radius, r, is then given by 
(ZLS99), 



Fino(r) 

where 



3h(r) 
16tt 2 
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h(r) 



(4/3) x 
-r- 2 )E(r 2 ) 
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(2r - r" 1 ) E (r~ 2 ) + 2 (r^ 1 - r) K (r^ 2 )] , 



r < 1, 
r > 1, 



(A2) 



where E and K are complete elliptic integrals. The luminosity of 
the disc due to reemission of photons incident on it is 



Lin C (d) =4tt / drrF inc (r). 



(A3) 



This can be divided into the contributions to the integral from the 
parts of the disc at d < 1 and d > 1, L[ nc = L^ c + Lf^ c , where 
the first term is nonzero only for d < 1. 

Then, the relative strength of Compton reflection can be iden- 
tified (ZLS99) with the ratio of the luminosity of the disc due to 
irradiation to the fraction of the sphere luminosity that is not inci- 
dent on the disc. However, we correct here for attenuation of the 
reflection from the parts of the disc with d < 1 because of scatter- 
ing by the hot electrons in the sphere. For the radial optical depth 
of the sphere r, we can then write 



R ■ 



(A4) 



1 — Line 

Another effect not included in the treatment of ZLS99 is the 
intrinsic dissipation in the disc. Far away from the center (so any 
effect of the inner boundary condition is negligible), the dissipated 
flux per unit area is tx r~ 3 , and we assume it for r > 1. On the 
other hand, the dissipation in the part of the disc inside the hot 
sphere, r < 1, is reduced due to the transfer of the power to the 
hot plasma. We assume here that dissipation either to be null or 
constant matching that of the outside disk, 



„ / \ Li n t f or 1, r < 1, 
F **W = -to\r->, r>l, 



(A5) 



where we parametrized the relative intrinsic dissipation by the di- 
mensionless factor, Lint, defined as the intrinsic luminosity of the 
disc extending from r = 1 to infinity (regardless of the actual value 
of d). 

The power in seed photons scattered in the sphere from both 
the reprocessing and the intrinsic dissipation (assuming r = 1) is 
then 



L s {d) 



3(1 -a) 
4tt 2 



dr rh (r) H — 

7T 



dr 



fc(r) 



max(d, 1) 



+ — ^ / drrh(r) 



(A6) 



where a is the albedo and the term given by the second line of 
equation <A6> appears only in the case of the nonzero dissipation 
within the sphere and when d < 1. 

The amplification factor of the process of thermal Comp- 
tonization is then A(d) = 1/L S (d) . It can be related to the spectral 
index of the power-law part of the Comptonization spectrum using 
e.g. the formula of MBP01 

r^C^-l)" 4 , (A7) 
where C = 2.19 and 8 = 0.14 for BHB with fcT sccd = 0.2 keV. 



